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In two dimensions strongly interacting bosons in a magnetic field can realize a bosonic integer 
quantum Hall state, the simplest two dimensional example of a symmetry protected topological 
phase. We propose a realistic implementation of this phase using an optical flux lattice. Through 
exact diagonalization calculations, we show that the system exhibits a clear bulk gap and the 
topological signature of the bosonic integer quantum Hall state. In particular, the calculation of the 
many-body Chern number leads to a quantized Hall conductance in agreement with the analytical 
predictions. We also study the stability of the phase with respect to some of the experimentally 
relevant parameters. 


Introduction — Symmetry Protected Topological 
(SPT) insulating phases of interacting bosons are the 
analogues of the celebrated free fermion topological in¬ 
sulators for bosonic systems (see Ref. [T| for a short re¬ 
view). The SPT phases have a bulk gap but no intrin¬ 
sic topological order. As a consequence, they do not 
host any fractionalized excitations. Still they can exhibit 
protected boundary excitations. Non-interacting boson 
phases are topologically trivial. The need for strong in¬ 
teractions to obtain SPT phases makes them harder to 
study than their fermionic cousins, where the band topol¬ 
ogy is the only required information. Even their coho¬ 
mology classification [2HB] does not cover all the known 
examples [7l412j . Therefore, finding experimentally rele¬ 
vant microscopic models in dimension larger than one 
remains an important and difficult task. 

The prototype of a SPT phase in two dimensions is the 
boson integer quantum Hall (bIQH) state. The physical 
properties of this state have been studied in Refs. IT51 - [TS1 
In particular, the Hall conductivity was shown [13] to be 
quantized and equal to an even integer. The edge physics 
consists of a charged chiral edge mode and a counter 
propagating neutral mode. Despite being non-chiral this 
edge structure is protected so long as charge conserva¬ 
tion symmetry is preserved. Following the proposal of 
Ref. HU several recent numerical studies [TTHTU] have a 
pointed out the possible emergence of bIQH in bilayer 
bosonic fractional Hall system at filling factor v = 1 for 
each component. In particular, a robust gapped phase 
has been observed, and some indications of the topolog¬ 
ical nature has been obtained. While providing a proof 
of concept, this setup would be difficult to realize in an 
experimental ultracold atomic system. Even the numeri¬ 
cal simulations still lack direct evidence of the topological 
nature of the observed state, such as the existence of edge 
modes or the measurement of the Hall conductance. 

Motivated by experimental settings in ultracold atomic 


gases, lattice versions with “Chern bands” could replace 
the lowest Landau level for the bIQH while providing 
a realistic implementation. In such systems, the emer¬ 
gence of strongly correlated topological phases has been 
studied in both the context of the fractional quantum 
Hall (FQH) effect in the presence of a lattice [20H25] or 
in a fractionally filled Chern band, namely a fractional 
Chern insulator [25H55] . As opposed to a single Landau 
level, lattice models can be engineered to have bands 
with a Chern number C higher than one. Since in the 
non-interacting regime these bands can be viewed as C 
copies of a Chern one band, a single C = 2 band would 
mimic a bilayer system. The fate of a partially filled 
band carrying a Chern number C > 1 in the strongly in¬ 
teracting regime has been recently studied [34MT] . Some 
of these systems host new phases that are generaliza¬ 
tions of the Halperin states [32] in the FQH with color- 
orbit couplings. It has also been shown that the C = 2 
band of the Harper-Hofstadter model can support a non- 
fractionalized bosonic phase[25], 4L at least at low par¬ 
ticle densities. Still, a simple lattice model of a Chern 
insulator in which the bIQH state can appear at high 
particle densities remains lacking. 

In this letter, we consider the implementation of the 
bIQH in optical flux lattices [43ti46] as a potential ex¬ 
perimental realization in ultracold atomic gases. Opti¬ 
cal flux lattices can be designed in several ways, offering 
a large control over the band topology, dispersion and 
Berry curvature distribution through parameters such as 
the number of internal states and the laser couplings |45j . 
These systems also allow to tune the Chern number of the 
lowest band while preserving its approximate flatness. In 
presence of strong interaction, we have shown the emer¬ 
gence of Halperin-like states^. Thus they constitute 
a natural candidate to look for the bIQH. Using exact 
diagonalization, we provide compelling evidence for the 
emergence of such a phase when fully filling the lowest 
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band of a lattice for which the Cliern number has been 
set to two. In particular, we compute the many-body 
Chern number, thus showing the quantization of the Hall 
conductivity. 

We consider bosonic atoms with N internal degrees of 
freedom in the optical flux lattice model introduced by 
Ref. 35] The one-body Hamiltonian is given by 

+ (1) 

where Ijy is the N x N identity matrix. The coupling 
of the different internal degrees of freedom using two- 
photon Raman transitions is described by the poten¬ 
tial V(r), with a characteristic scale V set by the laser 
strength. The laser beams are arranged to form a tri¬ 
angular lattice, and induce a set of allowed momentum 
transfers in reciprocal space with a triangular pattern. 

The characteristic kinetic energy scale in the Hamilto- 
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nian is the recoil energy E R = where k is the wave 
number of the laser beams. 

A synthetic gauge field in reciprocal space can be cre¬ 
ated by controlling the relative phases between the dif¬ 
ferent laser beams. These phases can be chosen so that 
the lowest band possesses any Chern number C smaller 
than N, the ratio C/N controlling the band dispersion. 
In the following, we will focus on C = 2. More details 
on the one-body model can be found in Ref. [J5J 0Bj. In 
Fig#, we give a schematic description of the lattice in 
reciprocal space. It is a triangular lattice spanned by the 
vectors m = (1, 0)/c and n 2 = yy k. Due to a higher 
translation symmetry the first Brillouin zone is spanned 
by Nki and k 2 . As a result, the real space unit cell 
spanned by a\ = -^(1, — 1/\/3) and a 2 = ^=(0,1), has 
an aspect ratio equal to N. A typical density of states is 
shown in Fig. [T]>. It clearly shows the flatness of the low¬ 
est band, quantified by its spread dib, and the large gap 
Aib to the second band. In this letter, we will consider 
the two numbers of internal degrees of freedom N = 5 
and N = 6. At the level of the one-body spectrum, we al¬ 
ready observe in Fig. [lj: a sharp difference between these 
two cases with respect to the laser strength. Indeed, the 
gap is rapidly decreasing with V for N = 5 while it stays 
mostly unchanged for N = 6. 

While Chern insulators are defined by a tight-binding 
model, optical flux lattice models are continuous in real 
space. Thus, the interaction between the atoms mainly 
depends on the considered atomic species. Here, we fo¬ 
cus on the simplest and more realistic interaction: the 
s-wave scattering that correctly describes cold gases of 
alkali atoms like 8 'Rb. Thus, the interaction potential is 
given by 

H in t = V int S(r - r') . (2) 

Previous studies [TTHTU] on the emergence of the bIQH 
within a bosonic bilayer FQH system involved interac- 
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FIG. 1. (a) Schematic description of the first Brillouin zone 
of the OFL model that we have considered for N = 5. The 
blue, red and black links represents the hopping amplitudes 
defined by the potential V (r). The phases are set to produce 
a lowest band carrying a Chern number C = 2. (b) Density of 
states for V/Er = 3 and N = 6. Note that we have rescaled 
the density of states of the lowest band by a factor 0.1. The 
band spread of the lowest band is <5ib = 0.0126 Er and the 
band gap between the lowest band and the second band is 
Aib = 0.721 Er. (c) Band spread and band gap of the lowest 
band as a function of the laser strength V for N = 5 and 
N = 6. 


tions both within each layer and between the two lay¬ 
ers, also given by the interaction potential of Eq. [2] A 
major difference in that case is the ability to tune the 
ratio of the interaction strengths between the layers and 
within a layer. Indeed, the bIQH was found[19] for ratios 
between 0.8 and 1.3. While a Chern C = 2 band can 
be decoupled into two copies of a Chern one band, this 
cannot be achieved once interactions are included. As a 
consequence, the notion of intra and inter layer interac¬ 
tion is not meaningful in this context and the interaction 
strength ratio is locked to one [551. 

We now turn to the numerical study of the model. We 
consider N R bosons on a finite size system with periodic 
boundary conditions defined by the vectors L x = N x ai 
and L y = N y a 2 . The number of orbitals per band is then 
equal to N x N y and, thus, the filling factor is v = jrfr- 
Here, we focus on v equal or close to 1, i.e. one boson 
per state of the lowest C = 2 band. 

Similar to the lowest Landau level projection in the nu¬ 
merical studies of FQH effect, we project this interaction 
onto the lowest band of the one-body model, removing 
the effect of both band mixing and band dispersion. Such 
an approximation is valid as long as <5ib -C Vj n t <C Aib- 
At unit filling of fermions, this approximation would lead 
to the usual integer quantum Hall (IQH) effect, the in¬ 
teraction only shifting the system energy. This is not 
case for the bosonic atoms considered where the non- 
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interacting and strong interacting regimes differ. 

A typical energy spectrum at unit filling is she 
Fig. [2jr. It exhibits a single ground state with a zero mo¬ 
mentum separated from the excited states by an energy 
gap A. This is expected for the bIQH : having no intrinsic 
topological order, its degeneracy should be equal to one 
irrespective of the genus of the surface we consider, in¬ 
cluding our torus geometry. Other phases that have been 
proposed in such a system, like two copies of the Moore- 
Read state [?7] or the non-abelian spin singlet state [38], 
would exhibit a non-trivial degeneracy. Note that we ob¬ 
serve this single state at zero momentum separated by 
an energy gap irrespective of the particle number, show¬ 
ing the absence parity effect: in a bilayer setup, Ng is 
required to be even by construction. 

Another feature that can be tested is the number of 
quasiholes states. For Chern insulators, this can be ac¬ 
cessed by increasing the number of sites compared to the 
ground state at fixed number of particles or by removing 
particles at fixed number of sites. Here, we choose the 
latter solution and show in Fig. ih the energy spectrum 
with one particle less than in Fig. UK This spectrum 
exhibits an almost flat lowest energy band consisting of 
one level per momentum sector. This is analogous to the 
IQHE where the number of states when adding a quasi¬ 
hole is equal to the number of orbitals in each Landau 
levels. 

To show that the features of the bIQH persist in this 
system in the thermodynamic limit, we show the evolu¬ 
tion of the many-body gap A as a function of system size 
in Fig. [3] for TV = 5 and N = 6. While the gap exhibits 
some fluctuations resulting from the aspect ratio varia¬ 
tion (a known effect[3j3l ;5Di for small sizes in any frac¬ 
tional Chern insulator), it shows a tendency to converge 
for larger systems around A ~ 0.1I^ nt . Note that we 
have also considered lower values of N such as N = 4. 
There size effects are even more pronounced, including 
the absence of a gap in some cases. 
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FIG. 2. Left Panel: Low energy spectrum of the two body 
interaction at v = 1 for N = 6, Nb = 14, N x = 7, N y = 3 
and V = 3Er. The energies are plotted as a function of the 
linearized momentum index where (K x ,K y ) are the two in¬ 
tegers defining the momentum sector. Note that the energies 
are shifted by the ground state energy Eq. Right Panel: En¬ 
ergy spectrum of the two body interaction at v = 1 for N = 6, 
Nb = 13, N x = 5, N y = 3 and V = 3 Er. Compared to the 
left panel, there is one particle less and thus one added quasi¬ 
hole. We observe an almost flat band of low energy states 
with one state per sector as expected for the IQH effect. 
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FIG. 3. Gap as a function of the inverse of the particle number 
for N = 5 (red dots) and N = 6 (blue triangles) and V = 3 Er 
up to Nb = 17. For each particle number, we choose N x and 
N y so that the aspect ratio is the closest possible to 1. How¬ 
ever, there is still important aspect ratio variations between 
the different point which prevent from doing an relevant ex¬ 
trapolation of the gap at the thermodynamical limit. Quite 
clearly, this gap does not close at the thermodynamical limit 
and will be of order 0.1 Vint 

To establish on a stronger footing the realization of the 
bIQH in this system, we compute the many-body Chern 
number of the ground state IT). A non-zero value of this 
number would show the quantization of the Hall conduc¬ 
tance, a direct signature of the bIQH. The many-body 
Chern number can be computed either via the contour 
integal 


C M b = ^Im(T|VT) 


(3) 


or, equivalently, using Stokes’ theorem, by the surface 
integral: 

(4) 

where 6 X and 8 y are the two fluxes that can be inserted 
through the two independent non-contractible loops on 
the torus. While these two integrals are equal, their nu¬ 
merical evaluation requires dividing the ( 0 x ,6 y ) torus in 
small patches of size 27r<5 x 2ttS. This discretization yields 
slightly different results. We show the many-body Chern 
number computed using both methods for a system of 
Ng = 12 bosons with N = 5 and N = 6 as a function 
of S in Fig. [4j Generically, the contour integral tends 
to overestimate the result whereas the surface one un¬ 
derestimates it. We have checked our numerical proce¬ 
dure against the half filled C = 1 band where Laughlin 
v = 1/2 physics is realized [3B] and we have found a 
similar behavior. These results are shown in the inset 
of Fig. [4j In both cases, we found that the Chern num¬ 
ber is equal to 2 within 1% precision. We have obtained 
similar results for all systems from Nb = 9 to Nb = 12, 
irrespective of the parity of Ng. This is an unambiguous 
evidence of the Hall conductance quantization and thus 
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the realization of the bosonic integer quantum Hall effect 
in this system. 



FIG. 4. Many-body Chern number Cmb computed both using 
the surface integral and the path integral version, as function 
of the discretization step 5 of the ( 9 x ,8 y ) torus for N = 5 
and N = 6. The system considered here has Ng = 12, N x = 
6,N y = 2 and V = 2>Eg. The dashed horizontal line at Cmb = 
2 is a guide for the eye. In the inset, we show the result at 
v = 1 when the lowest band has C = 1 where Laughlin state 
is realized for Ng = 6, N x = 4, N y = 3 and N = 5. 

Until now, we have focused on a specific set of param¬ 
eters for our model. It is experimentally relevant to look 
at the effect of the laser strength or the flat band approx¬ 
imation on the bIQH phase stability. We first consider 
the effect of the laser strength in the flat band approxi¬ 
mation. The numerical results for Nb = 14 are shown in 
Fig. [5] Remarkably, we observe that the many-body gap 
has the same behavior as the one-body gap that we 
have computed in Fig. [TJ:. The maximum for N = 6 is in 
the regime V — Er of a shallow optical lattice, empha¬ 
sising that this system is far from the tight-binding limit. 
A more important question is the fate of the bIQH phase 
away from the flat band approximation. Due to the nu¬ 
merical limitation, we still have to truncate the one-body 
basis. We thus consider the two lowest bands including 
their dispersion relation and look at the effect of the in¬ 
teraction strength Vint- The many-body gap is shown in 
Fig-! We have not directly computed the Chern num¬ 
ber in this case. However, as interaction strength V5nt in¬ 
creases and band mixing becomes more relevant, we find 
smooth evolution, with no gap collapse or phase transi¬ 
tion. Thus, these states are adiabatically connected to 
the states formed in the projected lowest band, so must 
retain the same Chern number C = 2 that we computed 
in that case. We have also checked that we can continu¬ 
ously deform the system up to the flat band limit without 
gap closing, as shown in the inset of Fig. [6j 

Conclusion — In this letter, we have considered the 
realization of the C = 2 bosonic integer quantum Hall 
effect in the lowest band of an optical flux lattice. We 
have characterized the topological order of the phase 
through the quashihole excitations and the evaluation 
of the electrical Hall conductance via many-body Chern 
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FIG. 5. Influence of the laser strength V on the many-body 
gap A for Ng = 14, N x = 7 and N y = 2. We have considered 
both N = 5 (red dots) and N = 6 (blue triangles). The 
behavior of A exactly mimics the behavior of the one-body 
gap Aib shown in Fig. [TJ:. 
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FIG. 6. Many-body gap as a function of the interaction 
strength Vint when the two lowest bands and their disper¬ 
sion relations are taken into account, for N = 6, Ng = 10, 
N x = 10, N v = 1 and V = 3 Eg. For comparison, the one- 
body gap is Aib = 0.721 Eg. The inset gives the continuous 
evolution at Vint = 0.175 -Er from the two band model (A = 0) 
to the flat band limit (A = 1). Both evolutions of the gap 
(black line) and the overlap with the flat band ground state 
(green line) show the adiabatic connection between these two 
systems. 


number. 

Note: During the final stage of this letter, we became 
aware of a related paper Ref. m providing numerical 
evidence for the bIQH in another lattice model with 
correlated hopping and a background gauge field. 
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